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ABSTRACT 

We present initial results from the first 3-diniensional numerical magnetohydrody- 
namical (MHD) simulations of magnetic field evolution in merging clusters of galaxies. 
Within the framework of idealized initial conditions similar to our previous work, we 
look at the gasdynamics and the magnetic field evolution during a major merger event 
in order to examine the suggestion that shocks and turbulence generated during a clus- 
ter/subcluster merger can produce magnetic field amplification and relativistic particle 



acceleration and, as such, may play a role in the formation and evolution of cluster- 
wide radio halos. The ICM, as represented by the equations of ideal MHD, is evolved 
self-consistently within a changing gravitational potential defined largely by the colli- 
sionless dark matter component represented by an N-body particle distribution. The 
MHD equations are solved by the Eulerian, finite-difference code, ZEUS. The particles 
are evolved by a standard particle-mesh (PM) code. We find significant evolution of 
the magnetic field structure and strength during two distinct epochs of the merger evo- 
lution. In the first, the field becomes quite filamentary as a result of stretching and 
compression caused by shocks and bulk flows during infall, but only minimal ampliflca- 
tion occurs. In the second, ampliflcation of the field occurs more rapidly, particularly 
in localized regions, as the bulk fiow is replaced by turbulent motions (i.e., eddies). 
The total magnetic field energy is seen to increase by nearly a factor of three over that 
seen in a non-merging cluster. In localized regions (associated with high vorticity), the 
magnetic energy can increase by a factor of 20 or more. A power spectrum analysis of 
the magnetic energy shows the amplification is largely confined to scales comparable 
to and smaller than the cluster cores indicating that the core dimensions define the 
injection scale. Although the cluster cores are numerically well-resolved, we cannot 
resolve the formation of eddies on scales smaller than approximately half a core radius. 
Consequently, the field amplification noted here likely represents a lower limit. We dis- 
cuss the effects of anomalous resistivity associated with the finite numerical resolution 
of our simulations on the observed field amplification. 

Subject headings: magnetohydrodynamics - methods: numerical- galaxies: intergalac- 
tic medium - galaxies: clusters: general 



1. INTRODUCTION 

Most radio sources in clusters of galaxies are 
associated with individual galaxies. There is, 
however, a class of radio sources known as radio 
halos (Jaffe 1977; Hanisch 1980, 1982) which ap- 
pears to be intrinsic to the cluster itself rather 
than any particular galaxy. Large-scale (>0.5 
Mpc) radio halos are rare, diffuse sources which, 
by their synchrotron emission, demonstrate the 
existence of magnetic fields (and a population of 
relativistic particles) on megaparsec scales. Here, 
we examine numerically the evolution of cluster- 
wide magnetic fields. 

It has been suggested that the magnetic field 
is generated by a turbulent dynamo mechanism 
(Jaffe, 1980; Roland 1981; Ruzmaikin, Sokoloff 
& Shukurov 1989; Burns et al. 1992; Tribble 
1993). Roland (1981) further suggested that 
galaxy wakes may be responsible for generating 
turbulence throughout the cluster. However, this 
model has difficulty explaining the radio observa- 
tions (Tribble 1991; Goldman & Raphaeli 1991; 
DeYoung 1992) because there does not appear 
to be enough energy in galaxy wakes to power 
the radio source. It has been noted that clusters 
containing radio halos (such as Coma, A2255, 
A2256, A2163 etc.) all contain evidence of recent 
dynamical evolution. That is, they all exhibit 
substructure such as non-Gaussian galaxy dis- 
tributions, multiple X-ray peaks, non-isothermal 
temperature distributions, and they lack cool- 
ing flows (Edge, Stewart, & Fabian 1992; Watts 
1992; Burns et al. 1995; Burns 1998). This 
suggests that cluster mergers are responsible, at 
least in part, for the formation of radio halos and, 
by inference, the growth of magnetic fields on 
cluster scales (De Young 1992; Tribble 1993). 

Cluster mergers are capable of supplying large 
amounts of kinetic energy, comparable to the 
thermal energy of the cluster, over megaparsec 
scales. Previous numerical simulations (Roet- 
tiger. Burns & Loken 1993; Schindler & Miiller 
1993; Pearce, Thomas & Couchman 1994; Roet- 



tiger, Loken & Burns 1997) show that merg- 
ers generate shocks, bulk flows and turbulence 
within the ICM. The first two of these processes 
can result in some field amplification simply 
through compression. However, it is the turbu- 
lence which is the most promising source of non- 
linear amplification. Provided the turbulence has 
non- vanishing helicity ( ( v-Vx v ) 7^ 0), weak 
fields can be amplified exponentially via the well- 
known a-effect (e.g., Ruzmaikin et al. 1989). 
Once the field becomes strong enough to have a 
significant back reaction on the flow, the system 
may enter an MHD dynamo regime. However, 
because the full system of MHD equations must 
be solved to study this regime, the nature and 
extent of field amplification which occurs in the 
MHD dynamo is not well understood. Of course, 
in order to account for the synchrotron emission, 
not only must there be field amplification, but 
also the merger induced magnetohydro dynamics 
must reaccelerate the relativistic particle popula- 
tion (Pacholczyk & Scott 1976; Eilek & Henrik- 
sen 1984). Whether the turbulence is responsible 
for both, or whether shocks alone are largely re- 
sponsible for particle acceleration is also not well 
understood. 

In this paper, we begin a study of magnetic 
field evolution in clusters of galaxies using three 
dimensional direct numerical MHD simulations. 
Here, we focus in particular on the role of clus- 
ter/subcluster mergers, the resulting gasdynam- 
ics (shocks, bulk flows, turbulence), and its role 
in the evolution (structure and amplification) of 
the cluster-wide magnetic field. Our numerical 
method is similar to the cluster merger studies 
that we have performed in the past, except in this 
case, we use an ideal MHD rather than a hydro- 
dynamics code to evolve the ICM. As in previous 
simulations, we use idealized, non-cosmological 
initial conditions consisting of two isothermal 
spheres in hydrostatic equilibrium which merge 
due to the influence of their mutual gravitational 
attraction. This allows us to study in detail the 
effects of the merger itself separate from the gen- 



eral cluster formation process thereby optimiz- 
ing numerical resolution through an efficient use 
of the computational volume. It also allows us 
to control both the structural and merger pa- 
rameters of the two clusters, thus allowing for a 
systematic survey of parameter space. Since the 
cluster magnetic field parameters (spatial distri- 
bution, strength, length-scale, etc.) are so poorly 
constrained at this time, we feel a parameter sur- 
vey is absolutely essential and will be the sub- 
ject of a future study. In this initial study, we 
present the results of only a single merger per- 
formed at two numerical resolutions, and a single 
non-merging cluster. 

In §|2|, we discuss our numerical method. Our 
initial conditions are presented in ^. A discus- 
sion of the magnetic field evolution, including nu- 
merical resolution effects, is presented in ^. We 
summarize our results in §^. 

2. NUMERICAL METHOD 

The ICM and magnetic fields therein are evolved 
using ZEUS (Stone & Norman 1992a,b), an Eu- 
lerian, finite-difference code which solves, self- 
consistently, the equations of ideal magnetohy- 
drodynamics (Jackson 1975). The numerical evo- 
lution of the magnetic field components is per- 
formed by the constrained transport (CT) algo- 
rithm (Evans & Hawley 1988) which guarantees 
preservation of the divergence-free constraint at 
all times. The method of characteristics (MOC) 
is used for computing the electromotive force 
(Hawley & Stone 1995). An extensive series of 
MHD test problems have demonstrated that the 
MOC-CT method provides for the accurate evo- 
lution of all modes of MHD wave families (Stone 
et al. 1992). We employ outflow boundary con- 
ditions on the MHD. 

The collisionless dark matter is evolved using 
an N-body code based on a standard particle- 
mesh algorithm (PM, Hockney &: Eastwood 1988). 
The particles and gas are evolved on the same 
grid using the same time step. The time step 



is determined by applying the Courant condi- 
tion simultaneously to both the dark matter and 
the magnetohydrodynamics. The only interac- 
tion between the collisionless particles and the 
gas is gravitational. Since we are modeling an 
isolated region, the boundary conditions for Pois- 
son's equation are determined by a multipole ex- 
pansion (Jackson 1975) of the total mass distri- 
bution (dark matter and gas) contained within 
the computational grid. Particles that leave the 
grid are lost to the simulation. Typically, less 
than a few percent of the particles leave the grid. 
The hybrid ZEUS/PM code was parallelized 
using the Message Passing Interface (MPI; Gropp, 
Lusk and Skjellum 1994). These simulations 
were run on the Cray T3E in the Earth and Space 
Data Computing Division of the NASA Goddard 
Space Flight Center. 

In this paper, we report the results of three 
simulations. In order to assess the effects of nu- 
merical resolution, we have conducted two simu- 
lations in which the resolution differed by a fac- 
tor of four. The low resolution simulation was 
performed on a uniform mesh with dimensions 
256 X 128^ zones. The high resolution simula- 
tion has dimensions which are a factor of two 
larger (512 x 256^ zones), but we gained an ad- 
ditional factor of two in resolution (over the cen- 
tral sub- volume occupied by the cluster cores and 
the merger axis) by implementing a non-uniform 
mesh. The high resolution mesh is uniform from 
zone 100 to 412 along the merger axis and from 
zone 90 to 166 perpendicular to the merger axis. 
Outside of this central region, we gradually in- 
crease the zone sizes, by ~3%, from one zone to 
the next radially outward. The physical dimen- 
sions of the computational volume are scaled to 
approximately 12 x 6.5 x 6.5 Mpc. For the phys- 
ical scaling describe in §0, our maximum resolu- 
tion corresponds to 50 kpc (4.3 zones/core ra- 
dius) and 12.5 kpc (17.2 zones/core radius), re- 
spectively, for the low and high resolution simu- 
lations. The third simulation, also conducted at 
high resolution, was designed to assess the nu- 



merical dissipation of the magnetic field in the 
absence of a merger. Here, we simply evolved a 
single cluster allowing it to move across the grid 
unperturbed. This simulation provides the base- 
line for comparison of magnetic field evolution in 
the merger. During the single cluster evolution, 
we note an initial rapid dissipation of the total 
magnetic field energy, ~50%, after which we find 
variations of less than 20%. 

3. INITIAL CONDITIONS 

3.1. Dark Matter and Gas 

Our initial conditions are similar to those used 
in our previous studies (e.g., Roettiger, Burns 
& Loken 1993; Roettiger, Loken & Burns 1997; 
Roettiger, Stone k Mushotzky 1997, 1998). We 
begin with two clusters whose gas distributions 
are consistent with observations of relaxed sys- 
tems. It can be argued that our dark matter 
profiles which, like the gas profiles, exhibit a flat 
core are inconsistent with recent strong lensing 
observations. These observations seem to im- 
ply a central cusp in the dark matter distribu- 
tion which may be more consistent with the dark 
matter distribution derived by Navarro, Frenk & 
White (1997) from cosmological simulations of 
structure formation. We do not believe that this 
difference will significantly alter the merger dy- 
namics. Both clusters in our simulations have 
mass distributions based on the lowered isother- 
mal King model described in Binney & Tremaine 
(1987). The lowered isothermal King model is a 
family of mass distributions characterized by the 
quantity ip/cr"^ which essentially defines the con- 
centration of matter. As ip/cr'^ increases, the core 
radius (re) decreases with respect to the tidal ra- 
dius (rt). We have chosen a model with ip/a'^=12 
in which we have truncated the density distribu- 
tion at 15rc- Near the half-mass radius, the to- 
tal mass density follows a power-law distribution, 
p r^ r^", where a ^^ 2.6. This model is consis- 
tent with mass distributions produced by cosmo- 
logical N-body simulations which show a ^ 2.4 



in a high density universe (^2=1) and ~2.9 in a 
low density universe ($7=0.2) (Crone, Evrard, & 
Richstone 1994). Initially the gas distribution 
(ICM) is in hydrostatic equilibrium within the 
gravitational potential defined by both the gas 
and dark matter components. The gas distri- 
bution is isothermal within the central drc- At 
larger radii, the temperature drops gradually. 

Our simulations are conducted in arbitrary 
units which can be scaled to physically mean- 
ingful parameters by chosing a mass and length 
scale. Table 1 contains the parameters of the 
two merging clusters scaled such that the more 
massive (or primary) cluster is representative of 
a rich Coma-like system. To date, radio halos 
have only been identified in relatively rich sys- 
tems. In line with the physical scaling used in 
Table 1, the density peaks are initially separated 
by 4.6 Mpc. Each cluster was given a small initial 
velocity in order to speed up the merger process 
and thus conserve computational resources. The 
final impact velocity (~2300 km s^"*^ ) is not af- 
fected significantly by the initial velocity which 
has components of 110 km s^^ parallel to the line 
of centers and 100 km s^^ perpendicular to the 
line of centers. These initial conditions result in 
a slightly off-axis merger with an impact param- 
eter less than 0.5rc. 

3.2. Magnetic Field 

Studies of cluster radio halos (Hanisch 1982) 
indicate field strengths of order 1 ^uG (e.g.. Coma, 
Kim et al. 1990; Ensslin & Biermann 1998; 
A2255, Burns et al. 1994; A2256, Rottgering et 
al. 1994) based on both Faraday rotation mea- 
sures and equipartition arguments. Unfortu- 
nately, neither of these methods is a direct mea- 
sure of the field strength, and they usually re- 
quire assumptions regarding several poorly con- 
strained parameters (e.g., Miley 1980). The only 
direct measurement of magnetic field strengths 
from synchrotron emission is via inverse Comp- 
ton scattering (Felton & Morrison 1966; Harris 
& Grindlay 1979). This method has been used 



by Bagchi, Pislar &; Lima Neto (1998) to con- 
strain the magnetic field strength in the diffuse 
steep spectrum radio source, 0038-096, in Abell 
85. They also find a mean field strength of ~l/iG. 
Other analyses of inverse Compton emission have 
indicated fields of 3-4/iG (Kaneda et al. 1995; 
Tashiro et al. 1996), but these appear to be as- 
sociated with individual radio galaxies and not 
the cluster itself. Radio sources in cooling flows 
have been shown to exhibit large Faraday rota- 
tion measures indicating field strengths in excess 
of 20/iG (Ge & Owen 1993; Owen & Eilek 1998), 
but, again, these observations are extremely local 
and may say more about the cooling flow environ- 
ment than about the cluster-wide fields (Soker & 
Sarazin 1994). 

The distribution of magnetic pressure on clus- 
ter scales is similarly uncertain. Deiss et al. (1997) 
show that the synchrotron halo in Coma, after 
subtracting point sources, traces the X-ray sur- 
face brightness distribution indicating that the 
magnetic pressure gradient is similar to that of 
the thermal pressure. In addition to the spatial 
distribution of the magnetic pressure, we must 
also be concerned with the spectral power dis- 
tribution. The magnetic fields are believed to 
exist as tangled flux ropes on a variety of scales 
(e.g., Ruzmaikin et al. 1989). The physical scale 
on which they are tangled is uncertain. Observa- 
tions of polarization in cluster radio sources indi- 
cate tangling on scales ranging from a kiloparsec 
(e.g., Feretti & Giovannini 1997) to tens or even 
hundreds of kiloparsecs (Kim et al. 1990; Kim, 
Tribble, & Kronberg 1991; Feretti et al. 1995) . 

Guided as closely by the observational data as 
possible, we begin by defining a magnetic vector 
potential, A(k)=Aok~'^, where the amplitudes 
(Ao) of each Cartesian coordinate are drawn from 
a Gaussian distribution. The vector potential 
is then transformed via a 3-dimensional Fast 
Fourier Transform (FFT) into physical space 
where it is scaled spatially by the gas density dis- 
tribution. Assuming a uniform spherical collapse 
(likely a gross over-simplification, e.g., Evrard 



1990; Bryan et al. 1994; among many others) 
and flux freezing, it can be shown that the mag- 
netic field will scale as pjas- Next, we scale the 
amplitude of the magnetic pressure such that the 
maximum magnetic energy density is equal to 1% 
of the local thermal pressure. This leads to a 
mean field of 0.07 //G within the central Ir^ of 
the primary cluster. Finally, we initialize a tan- 
gled divergence-free magnetic field from A via 
B = V X A. This results in B^ ex k^^^""). Here, 
we adopt a = 5/3. The premerger profile shapes 
(normalized at the core radius) can be seen in 
Figure |l|. 

4. DISCUSSION 

4.1. Magnetic Field Evolution 

A more detailed description of the hydrody- 
namical and N-body evolution of merging clus- 
ters of galaxies can be found in Roettiger, Burns, 
& Loken (1996); Roettiger, Loken, &; Burns 
(1997); Roettiger, Stone k Mushotzky (1997, 
1998). In fact, since the initial fields are weak, 
the hydrodynamics is largely the same. Here, we 
focus on the magnetic field evolution. 

Figure ^ shows the evolution of gas density 
(Column 1), gas temperature (Column 2) and 
magnetic pressure (B^, Column 3) at four epochs 
(Rows 1-4) during the merger. The epochs de- 
picted here correspond to 0.0, 1.3, 3.5, and 5.0 
Gyrs after the time of closest approach. Each 
panel represents a 2-dimensional slice taken through 
the core of the merger remnant and in the plane 
of the merger. The slice dimensions (3.75 x 3.75 
Mpc) are only a small fraction of the simulation 
volume. 

In Figure |2|a,b,c, the smaller of the two clus- 
ters (hereafter the subcluster) is seen impinging 
on the core of the larger cluster (or primary) from 
the right hand side and moving toward the left. 
Since the cluster gravitational potential is quite 
steep and since gas is continually being stripped 
from the subcluster, a gradient in the gas veloc- 
ity appears across the subcluster. The result is a 



stretching, and consequently, a small amplifica- 
tion of the magnetic field which produces long 
radial filaments, particularly in the subcluster 
wake (Figure ^). As the merging cores become 
coincident, the gravitational potential reaches an 
extreme minimum drawing in gas from all direc- 
tions which further enhances the radial filamen- 
tary structure. A shock forms along the leading 
edge of the subcluster resulting in compression 
and amplification of the magnetic field along the 
shock front, which is evident by the arc of hot 
gas visible in the temperature data. Figure |2|b. 

At 1.3 Gyrs (Figure ^,e,f), the bow shock has 
propagated off the left hand side of the frame 
having left behind a sheet of enhanced magnetic 
field which had been compressed along its lead- 
ing edge (Figure |2|f). In the subcluster wake, the 
magnetic field is drawn out of the primary cluster 
into long filaments. It is this type of magnetic 
field distribution which may explain the south- 
west extension of the Coma cluster radio halo 
(Deiss et al. 1997). If so, this would tend to 
support the assertion that the NGC 4839 group 
has already fallen through the core of Coma and 
is dragging magnetic field with it (Burns et 
al. 1994). 

Once the dark matter component of the sub- 
cluster exits the primary core, the gravitational 
potential minimum quickly returns to near pre- 
merger values. This causes a rapid expulsion of 
the gas that was drawn in during the merger. 
The expelled gas interacts with residual infalling 
gas from the subcluster creating a second shock 
that propagates upstream (to the right) result- 
ing in further compression and amplification of 
the field. Morphologically, there is considerable 
similarity between the two shocks and magnetic 
field structures (Figure ||f ) discussed here and the 
diffuse radio halo observed in A3667 (Rottgering 
et al. 1997). We address the possible physical 
connection between the merger induced magne- 
tohydrodynamics seen here and the A3667 radio 
halo in a future paper (Roettiger, Burns, &: Stone 
1999). 



During the early stages of the merger, the clus- 
ter gasdynamics are dominated by large scale 
bulk fiows and large eddies in the wake of the 
subcluster. Between 1.3 Gyrs and 2.5 Gyrs, the 
subcluster dark matter remnant passes through 
the cluster core three times. Although each pas- 
sage is less extreme than the one before, each 
contributes to the breakdown of the bulk flows 
and to the randomization of the gas velocities. It 
is during this phase that the most extreme am- 
plification of local field energy occurs (26 x over 
a non-merging cluster). After 2.5 Gyrs, the mag- 
netic field extrema begin to decay with time while 
the total magnetic energy continues to grow. The 
breakdown of the long filamentary structures is 
seen to begin in the core (Figure §i) and proceed 
radially outward as the cluster relaxes (Figure 

Figure y (solid line) shows the increase in total 
magnetic energy within the entire computational 
volume relative to the non- merging single cluster. 
Note the slight increase (~30%) during the first 
core passage (t =0). This is a result of the initial 
compression which, although it approaches the 
strong shock conditions, is very local and con- 
sequently does not greatly effect the total mag- 
netic energy within the larger volume. From core 
passage to the end of the merger simulation, the 
magnetic energy grows essentially linearly with 
time increasing by nearly a factor of 2.75 after 
4 Gyrs. We find that the fractional field ampli- 
fication within the cluster core is comparable to 
that seen in the much larger volume. Of course 
the initial compression at t =0 makes a much 
larger relative contribution to the smaller vol- 
ume. There is then a decline in the magnetic en- 
ergy as the field is drawn out of the cube only to 
return with the subcluster remnant after about 
1.3 Gyrs. The dashed line in Figure y shows the 
growth of the total magnetic field energy relative 
to the thermal energy also normalized to the sin- 
gle cluster evolution. After passage of the initial 
shock (i =0), the relative increase of magnetic 
energy to thermal energy tracks the relative in- 



crease in the magnetic energy quite well (i.e., the 
solid and dashed lines in Figure |3| are essentially 
parallel after t =0). Figure ^ shows the evolu- 
tion of the total thermal, kinetic and magnetic 
energy densities within 400 kpc of the gravita- 
tional potential minimum. These quantities are 
not normalized and have been scaled to the phys- 
ical dimensions in Table 1. Here we note the 
fractionally large increase in both thermal and 
kinetic energy, particularly at the time of core 
passage {t =0). The increase in magnetic energy 
is far more modest, indicating that a very small 
fraction of the merger energy goes into field am- 
plification. 

It is important to note that although the total 
magnetic energy has increased by only a factor 
of ~2.75, local amplification can be considerably 
greater. At ~2 Gyrs after the initial core passage, 
the maximum magnetic pressure has increased 
by a factor of 12 over premerger values and by a 
factor of 26 over the non-merging cluster at the 
same epoch. The peak maxima in magnetic pres- 
sure occur between 1.3 and 2.5 Gyrs. As men- 
tioned above, at 1.3 Gyrs, the initial bulk flows 
begin to breakdown with the second passage of 
the subcluster remnant, and at 2.5 Gyrs, stirring 
by the subcluster remnant largely ceases. The 
peak magnetic pressure after 5 Gyrs is still a fac- 
tor of three greater than both premerger values 
and the non-merging cluster at the same epoch. 

Figure |5| shows the azimuthally averaged and 
normalized profiles of gas density, thermal pres- 
sure and magnetic pressure at 5 Gyr after closest 
approach. Comparison with the initial profiles 
in Fig. H, shows relatively little evolution in the 
profile shapes. At this late stage of the merger 
the cluster properties have had time to equili- 
brate, and profiles tend to smooth out localized 
substructure. Numerical reconnection and dissi- 
pation of the magnetic field has proceeded in the 
core on scales less than 4-5 zones (50-60 kpc). 
We also note an increase in the thermal pressure 
of the core relative to the outer regions of the 
cluster. This is a common feature of the merger 



simulations (Roettiger et al. 1997). Much of the 
kinetic energy of the merger is dissipated within 
the core resulting in a rise in temperature accom- 
panied by an expansion of the gas distribution. 

Another way to look at the magnetic field 
evolution is through an analysis of the spectral 
power distribution. Figure ^ shows the ratio 
of the pre- to post-merger magnetic field power 
spectrum at three post-merger epochs. At 1.3 
Gyrs, before turbulent gas motions have replaced 
the bulk flows, power has increased by less than 
a factor of 10 on resolved (>4-5 zones) scales. At 
later times power has increased by nearly a fac- 
tor of 20 on the scale of ~4 zones and greater 
than a factor of 10 on scales less than 6-7 zones. 
Again, this is consistent with our factor of three 
increase in total magnetic energy since relatively 
little power resided in these small scales initially. 

Field amplification occurs on scales less than 
~8-10 zones. This is comparable to the diameter 
of the subcluster's gas core. Although numerical 
resolution is likely still playing a role, it is rea- 
sonable to expect that the largest scale on which 
turbulent eddies form and have sufficient time to 
turn over and amplify the field would be com- 
parable to the subcluster core dimensions. The 
important dimension is the gas core diameter, 
because it is the hydro dynamical interaction be- 
tween the clusters that is ultimately responsible 
for the field amplification, and the gas core sup- 
plies the greatest hydrodynamical impact. Sub- 
cluster gas at larger radii is much less dynami- 
cally significant (i.e., of much lower density and 
velocity) and is largely stripped by the time the 
cluster cores interact. Unfortunately, even at 8- 
10 zones we may still be uncomfortably close to 
the resolution limitations. Given greater resolu- 
tion, we would expect the eventual formation of 
smaller eddies which could potentially result in 
even greater field amplification on small scales. 
Small eddies turn over more rapidly than large 
eddies and will thus amplify the field to a greater 
extent over a given time scale. 



4.2. Resolution Study 

In order to study the effects of numerical res- 
olution, we have conducted two simulations that 
differed in effective resolution by a factor of four. 
In the low resolution simulation, the primary 
cluster core is resolved by 8-9 zones while in the 
high resolution simulation (the results of which 



were discussed in §4.1), the core is resolved by 
34-35 zones. In ZEUS, an artificial viscosity is 
used to thermalize kinetic energy in shocks; this 
viscosity smooths shocks over 4-5 zones. Simi- 
larly, the magnetic field resolution, (i.e., the scale 
over which numerical reconnection occurs) is ~4 
zones. In both simulations, we see the formation 
of long filaments during the early stages of the 
merger. The width of these filaments appear to 
be determined by the effective resolution. 

In the low resolution simulation, the field evo- 
lution largely ceases after the formation of the 
filaments. The filaments are advected about the 
cluster by residual gas motions, but no further 
amplification occurs. This is in contrast to the 
high resolution simulation which shows a factor 
of 2.75 increase in the total magnetic energy af- 
ter 4 Gyrs. It appears that the lack of numer- 
ical resolution prevents the development of tur- 
bulent eddies on meaningful scales. That is, the 
injection scale for turbulence in these mergers is 
comparable to the cluster core dimensions. Since 
the subcluster core is only marginally resolved 
in the low resolution simulation, turbulent ed- 
dies are not resolved well enough to form nor if 
they form do they persist long enough for ampli- 
fication of the magnetic field to occur. Figure ^ 
shows the ratio of the pre- to post-merger power 
spectra for both resolutions at 3.4 Gyrs after core 
passage. The low resolution simulation shows no 
significant amplification on scales greater than 4- 
5 zones, and the amplification that is evident is 
considerably less than it is in the high resolution 
simulation. At this time, there is no convergence 
in the magnetic field evolution. We speculate 
that given higher resolution we could follow the 
cascade of eddies to smaller and smaller scales 



expecting greater field amplification as the turn 
over timescale for the eddies decreases. Unfor- 
tunately, at this time, a significant increase in 
resolution over our high resolution simulation is 
not technically feasible given the limitations of 
currently available computational platforms. 

Numerical dissipation will have two effects on 
the MHD of the mergers: 1) numerical viscos- 
ity will artificially truncate the spectrum of tur- 
bulent eddies on scales of a few grid zones, and 
2) numerical resistivity will allow for anomalous 
dissipation of the field on similar scales. The ac- 
tual conductivity and viscosity of the ICM may 
have significantly different characteristic scales, 
so that the magnetic Prandtl number (the ra- 
tio of the coefficient of viscosity and resistivity) 
may be much different than one. Thus, study- 
ing the degree of amplification of the magnetic 
field as the magnetic Prandtl number is varied 
in the simulations is of great interest. However, 
since this requires a large dynamic range from 
core radius to the dissipation scales (the latter 
of which must be resolved in order to capture 
the turbulent eddies), such studies appear to be 
well beyond the range of current computational 
resources. Our largest simulation is far from hav- 
ing the requisite dynamic range in the remnant 
core for such studies. 

5. SUMMARY 

We have presented the initial results from 
3-dimensional numerical MHD/N-body simula- 
tions of merging clusters of galaxies. We find that 
cluster mergers can dramatically alter the local 
strength and structure of cluster-wide magnetic 
fields. Early in the merger, there is a filamen- 
tation of the field caused by tidal forces acting 
on the ICM component of both clusters. The 
steepness of the gravitational potential causes a 
gradient in the gas velocity which results in the 
stretching and filamentation of the field. We also 
find compression and amplification of the field 
along a shock located at the leading edge of the 
subcluster. Once formed, this structure prop- 



agates with the shock through the core of the 
cluster and out the other side. A similar fea- 
ture is formed on the upstream side of the merger 
when gas that is expelled from the rapidly vary- 
ing gravitational potential interacts with resid- 
ual infall from the subcluster. The result is two 
shocks, which through compression form mag- 
netic field structures on opposing sides of the 
merger remnant core. The strong shocks asso- 
ciated with these structures may be the source 
of energetic particles to power the synchrotron 
emission. These structures are similar in mor- 
phology to the radio halo in A3667 (Rottgering 
et al. 1997). Similarly, we find an extension of 
the magnetic field distribution along the merger 
axis toward the subcluster which is reminiscent 
of the southwest extension of the radio halo in 
Coma (Deiss et al. 1997). 

In the early stages of the merger, amplifica- 
tion of the field is limited to that produced by 
compression in the bow shock and by stretching 
resulting from infall. It is only after the merger 
induced bulk flows breakdown into turbulent gas 
motions that significant local field amplification 
takes place. Since we find basically two epochs 
of magnetic field evolution, we suggest that there 
are two epochs of radio halo formation. The first, 
represented by the halos in A2255 (Burns et 
al. 1995) and A2256 (Rottgering et al. 1994), are 
in the early stages of a merger and result largely 
from field compression associated with the ini- 
tial infall. The second, represented by the halo 
in Coma (Kim et al. 1990), is in a later stage 
of the merger dominated more by turbulent gas 
motions. 

We find that as a result of the merger the 
total magnetic energy within the computational 
volume rises steadily by nearly a factor of three 
during the 5 Gyrs after core passage. On smaller 
scales we find that the magnetic field energy in- 
creases by greater than a factor of 10 to 20. Of 
course, we have modeled only a single merger. It 
is likely that massive clusters will undergo sev- 
eral major mergers during their life time (every 



2-4 Gyrs; Edge et al. 1992) and that each succes- 
sive merger will further amplify the field. Also, 
it is likely that the ICM in massive clusters is be- 
ing stirred almost continuously by galaxy wakes 
(Roland 1981; De Young 1992; Merrifield 1998) 
and by residual infall including many lesser merg- 
ers (Norman & Bryan 1998). The increase in the 
magnetic energy is only about 10~^ of the ki- 
netic energy imparted by the merger event. The 
remaining energy goes into residual gas motions 
and heating of the ICM. This value may be lim- 
ited by numerical resolution effects. 

Our comparison of the premerger and post- 
merger power spectra of the magnetic field dis- 
tribution show that the greatest fractional ampli- 
fication of the field is on relatively small scales. 
There are several reasons for this. First, the 
most extreme amplification is associated with lo- 
calized gasdynamics such as shocks and small 
scale resolved eddies which have the fastest turn- 
over time and therefore the shortest amplifica- 
tion time. Second, the scale of the merger in- 
duced turbulence is strongest on scales smaller 
than the perturbation, in this case, the core di- 
ameter of the impinging subcluster (17 zones in 
the high resolution simulation). 

Resolution effects are always important in a 
study of this type. Observations indicate that 
magnetic fields, like turbulence, have power on a 
range of physical scales. It is impossible, at this 
time, to simulate all scales accurately. We note a 
couple of significant resolution effects. First, the 
width of the filaments formed during the early 
stages of the merger appears to be determined 
by the numerical resolution. This is also true 
of the magnetic sheets that form along shock 
structures. Second, and most importantly, in 
order to see significant amplification, turbulent 
eddies must be resolved. In these simulations 
eddies are resolved at 4-8 zones across. Finally, 
if the small scales (1 kpc; Feretti et al. 1995) 
for field tangling deduced from the radio obser- 
vations are correct, then we could be missing a 
very important dynamical scale in these simu- 
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lations. The dynamic range required to resolve 
the fuh spectrum of turbulent eddies on scales 
below the core radius, and yet still model the 
global merger, presents an enormous challenge 
to future numerical studies. Future work will in- 
clude a survey of magnetic field initial conditions 
(e.g., length scales, pressure distribution), and 
the modeling of specific clusters containing radio 
halos (e.g., A3667; Roettiger et al. 1999). 
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FIGURE CAPTIONS 

Fig. [l] Premerger gas density (solid line), 
thermal pressure (dotted line) and magnetic pres- 
sure (dashed line) profiles of the primary clus- 
ter normalized to their values at the core radius. 
The central dip in the magnetic pressure is the 
result of numerical reconnection of the field. At 
this time, the clusters have been evolved approx- 
imately 3 Gyrs. The merger will occur in approx- 
imately 1.5 Gyrs. 

Fig. ^ The evolution of log gas density (col- 
umn 1), gas temperature (column 2), and log 
magnetic pressure (column 3) in two-dimensional 
slices taken through the cluster core in the plane 
of the merger. Each row represents a different 
epoch. Row 1 is t=0, the time of core coinci- 
dence. Rows 2, 3, and 4 represent i=1.3, 3.4, 
and 5.0 Gyrs, respectively. In each quantity, 
black represents low values and red represents 
high values. The color transfer function is con- 
sistent from one panel to the next within a given 
physical quantity. Each panel is 3.75x3.75 Mpc. 
The range of values in each quantity are: gas 
density, 6xl0~^ to 9.2x10"^ cm~^; gas temper- 
ature, 2.6 to 23.1 keV; and magnetic pressure, 
1.24x10-21 to 5.9x10-13 dynes cm-^. 

Fig. ^ The total magnetic energy within the 
entire computational volume (solid line) and the 
evolution of the magnetic to thermal energy ratio 
(within the entire volume), also normalized to a 
non- merging cluster (dashed line). Time is given 
relative to the time of closest approach, or core 
passage. 

Fig. Q The total thermal (open diamonds), 
kinetic (open triangles) and magnetic (open squares) 
energy densities within an 800 kpc cube centered 
on the gravitational potential minimum. The 
quantities are scaled according to the physical 
dimensions in Table 1. Time is given relative to 
the time of closest approach, or core passage. 

Fig. ^ The gas density (solid line), ther- 
mal pressure (dotted line) and magnetic pressure 
(dashed line) profiles of the merger remnant nor- 



malized to their values at the core radius. The 
central dip in the magnetic pressure is the re- 
sult of numerical reconnection of the field. At 
this time, the clusters have been evolved approx- 
imately 5 Gyrs since core passage. 

Fig. ^ The ratio of pre- to post-merger 
magnetic energy power spectra at three differ- 
ent epochs, 1.3 Gyr (solid line), 3.4 Gyrs (dashed 
line), and 5.0 Gyrs (dash-dot line), for the high 
resolution merger. This analysis was performed 
on an 800 kpc cube (64'^) centered on the gravi- 
tational potential minimum. The scale for signif- 
icant amplification appears to be less than 8 to 
10 zones. 

Fig. The ratio of pre- to post-merger mag- 
netic energy power spectra for low (50 kpc/zone; 
solid line) and high (12.5 kpc/zone; dashed line) 
resolution simulations at 3.4 Gyrs. The analy- 
sis was performed on an 800 kpc cube centered 
on the gravitational potential minimum. Field 
amplification is significantly less in the low reso- 
lution simulation, and there is no significant am- 
plification on resolved scales (>4-5 zones) in the 
low resolution simulation. 
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Table 1. Initial Cluster Parameters 






Cluster 
ID 


ML 
(10^4 Mo) 


(keV) 


(km s-i) 


rt ft /5' 
(kpc) 


(10-3cm-3) 


„,8 
^impact 

(km s'l) 


Primary 
Subcluster 


8.0 
3.2 


6.7 
3.1 


785 
526 


220 0.12 0.75 
135 0.06 0.72 


1.55 
1.01 


2300 



Temperature. ^ Ve- 



1 Total Mass R<3 Mpc. 

locity Dispersion. 

^ Core Radius. ^ Global Gas Fraction, by mass. 

6 /3 = uruha^/kT; R<1.5 Mpc. 



^ Central Gas Density, 
matter). 



Impact Velocity (dark 
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